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基础 研究 
基于 动力 学 聚 类 与 a 散 度 测度 的 动态 心肌 PET 图像 因子 分 析 
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摘要 :目的 建立 一 种 基于 动力 学 聚 类 与 a 散 度 测度 的 因子 分 析 方 法 ,分 析 从 动态 心肌 PET 图 像 中 无 创 地 提取 血 输入 函数 及 局 部 
组 织 的 时 间 活 度 曲线 。 方 法 通过 最 小 化 动态 图 像 与 因子 模型 间 的 a 散 度 将 动态 PET 图 像 做 初步 的 分 解 , 得 到 初始 因子 与 因子 
图 像 ,然后 联合 PET 像 素 动力 学 聚 类 的 先 验 信息 解决 因子 分 析 模 型 中 解 的 非 唯一 性 问题 ,将 初始 因子 与 因子 图 像 通过 空间 变换 
生成 具有 生理 意义 的 组 织 活 度 曲线 及 组 织 空间 分 布 。 结 果 与 传统 的 最 小 二 乘法 测度 和 最 小 化 因子 图 像 间 重 苹 程度 约束 模型 
相 比 ,本 模型 对 噪声 的 敏感 性 较 低 ,提取 出 的 结果 的 精确 性 较 高 。 结 论 通过 选取 最 优 的 a 值 作为 因子 分 析 模 型 的 测度 ,并 引入 
PET 图 像 像 素 的 动力 学 聚 类 信息 ,能 精确 地 获得 血 输入 函数 及 局 部 组 织 的 时 间 活 度 曲 线 , 在 视觉 评价 及 量化 评价 均 具 有 优 
质 表 现 。 
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Kinetic cluster and a-divergence-based dynamic myocardial factorial analysis of positron- 


emission computed tomography images 
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Abstract: Objective We purpose a novel factor analysis method based on kinetic cluster and a-divergence measure for 
extracting the blood input function and the time-activity curve of the regional tissue from dynamic myocardial positron 
emission computed tomography(PET) images. Methods Dynamic PET images were decomposed into initial factors and factor 
images by minimizing the a-divergence between the factor model and actual image data. The kinetic clustering as a priori 
constraint was then incorporated into the model to solve the nonuniqueness problem, and the tissue time-activity curves and 
the tissue space distributions with physiological significance were generated. Results The model was applied to the "RbPET 
myocardial perfusion simulation data and compared with the traditional model-based least squares measure and the minimal 
spatial overlap constraint. The experimental results showed that the proposed model performed better than the traditional 
model in terms of both accuracy and sensitivity. Conclusion This method can select the optimal measure by a value, and 
incorporate the prior information of the kinetic clustering of PET image pixels to obtain the accurate time-activity curves of the 
tissue, which has shown good performance in visual evaluation and quantitative evaluation. 
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无 创 测定 局 部 心肌 血 流量 (RMBF) 一 直 是 临床 心 。 析 , 目 前 广泛 采用 的 估计 血 输入 函数 方法 为 基于 感 兴 
脏 学 所 追求 的 目标 之 一 ”。 通 过 从 动态 图 像 中 估计 血 区域 法 ,该 方法 简单 易 行 ,但 极度 依赖 于 医生 手工 勾画 
输入 函数 与 输出 函数 ,并 应 用 动力 学 模型 ,可 以 得 到 ” ”的 感 兴趣 区 域 ,而 感 兴 趣 区 域 的 准确 性 又 受 医生 个 人 经 
绝对 的 定量 参数 (如 物质 转运 速率 代谢 速度 和 受 体 验 及 部 分 容积 效应 影响 。 因 子 分 析 作 为 从 变量 群 中 提 
结合 率 等 功能 参数 等 ), 从 而 实现 对 心肌 血 流 的 定量 分 。 取 公 共 因 子 的 统计 技术 ”, 已 经 被 用 来 从 动态 图 像 中 无 
析 “™…。 在 动力 学 模型 的 应 用 中 , 血 输入 函数 的 估计 是 非 ” 创 地 提取 血 输入 函数 和 组 织 的 时 间 活 度 曲 线 ”", 如 基 
常 关键 的 ,结果 的 准确 与 否 直 接 影响 到 后 期 的 定量 分 。 于 主 成 分 分 析 ”*、 独 立成 分 分 析 ”“、 最 小 二 乘法 “"、 
最 大 似 然 法 “的 因子 分 析 模型 相继 被 提出 。 然 而 ,用 
子 分 析 方 法 获得 的 结果 存在 非 唯一 性 ,各 种 基于 先 验 信 
息 的 方法 被 提出 用 来 解决 结果 的 非 唯 一 性 “2 ,另外 ， 
Sitek 等 "提出 基于 最 小 化 各 因子 图 像 间 的 重叠 程 度 的 
优化 方法 可 以 获得 唯一 的 结果 。 

合理 地 构建 因子 模型 与 实际 数据 之 间 的 测度 是 ; 
确 建 模 的 关键 ,目前 在 因子 分 析 模型 中 常见 的 是 基于 
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斯 噪声 模型 的 最 小 二 乘 测度 与 基于 泊 松 噪声 模型 的 最 
大 似 然 测度 “*”。 而 临床 获取 动态 PET 图像 噪声 很 难 
直接 用 高 斯 或 泊 松 模型 进行 刻画 ,因此 ,本 文 提 出 基于 a 
散 度 的 因子 分 模型 ,a 散 度 由 S.Amari 提 出 ,是 诸多 信息 
散 度 族 中 的 一 种 ,通过 选取 不 同 的 a 值 ,可 以 用 于 任意 噪 
声 模型 的 数据 测度 2 。 

针对 因子 分 析 模型 中 解 的 非 唯一 性 ,传统 上 利用 最 
小 化 因子 图 像 间 的 重 秋 程度 这 一 约束 来 解决 该 问题 ,但 
由 于 受到 PET 图 像 分 状 率 和 部 分 容积 效应 的 影响 ,心肌 
信号 经 常会 摊 杂 着 10%~15% 的 血液 信号 ,基于 最 小 化 
子 图 像 之 间 的 重合 程度 的 约束 未 能 考虑 到 这 一 事 
实 。 在 动态 PET 人 研究 中 ,每 一 个 像素 的 时 间 行 为 可 以 通 
过 曲线 时 间 活 力 (TAC) 来 描述 ,我 们 将 具有 相同 时 间 行 
为 的 像素 进行 聚 类 ”” ,把 这 一 聚 类 结果 作为 先 验 约 束 用 
来 解决 因子 分 析 模 型 中 的 非 唯一 性 ,得 到 具有 生理 意义 
的 组 织 活 度 曲 线 。 因 子 分 析 模型 应 用 于 基于 XCAT 
(Extended Cardiac- Torso phantom) (Ki ^^ [Jj Cf "Rb 
PET 心 肌 灌注 仿真 数据 中 ,实验 均 表 明 我 们 的 方法 在 视 
党 评价 及 量化 评价 等 方面 优 于 传统 的 因子 分 析 模 型 。 


1 方法 
1.1 PET 图 像 的 获取 

对 于 动态 序列 PET 图 像 志 像素 点 /在 第 六 帧 的 示 踪 
剂 活 度 zw 可 通过 以 下 积分 式 计 算得 到 交 ， 


Dynamic sequence 


1 对 动态 序列 做 因子 分 析 示 意图 


Factor images 


Xin = fe. t)e"dt (1) 
其 中 ,c0,0 表 示 像 素 点 在 时 间 点 ! 的 示 踪 剂 活 度 ,如 
和 如. 分 别 表 示 第 症 帧 的 起 始 时 间 和 终止 时 间 。 那 么 
Xu Xp cto 
jx cx (2) 
x; X, Hind Xn, 
其 中 ren" 7 ,7 的 一 列表 示 一 帧 图 像 ,7 的 一 行 表 
示 单 个 像素 点 的 时 间 活 度 曲线 。 


动态 PET 序 列 图 像 可 分 解 为 有 限 的 时 间 序 列 向 量 
FE 和 与 其 相对 应 体 素 的 空间 分 布 C, 这 一 分 解 模型 可 用 
下 式 来 表示 

I=CF+e (3) 

其 中 动态 PET 图 像 1 的 大 小 为 Nx7, 其 中 入 表示 每 
一 帧 PET 图 像 的 像素 点 数目 ,7 表示 动态 网 像 的 帧 数 。 
因子 矩阵 Ff 的 大 小 为 Kx7, 因子 图 像 C 的 大 小 为 NxK,， 
其 中 天 为 所 需 提取 的 因子 数目 。s 表 示 剩 余 信 号 , 即 图 像 
的 噪声 和 被 模型 所 忽略 的 小 部 分 信号 。 在 这 里 ,为 保证 模 
型 的 实际 应 用 价值 ,因子 数目 K 要 满足 K<<N,K<<7。 

图 1 表示 对 动态 序列 PET 图 像 做 因子 分 析 的 示意 
图 ,在 理想 状况 下 ,因子 矩阵 正 对 应 于 动态 图 像 中 各 个 
生理 组 织 的 时 间 活 度 曲线 ,因子 图 像 C 对 应 于 各 个 生理 
组 织 的 空间 分 布 。 
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Fig.1 A schematic illustration of factor analysis for dynamic sequences. 


1.2 提出 的 因子 分 析 模 型 

子 分 析 模 型 的 求解 主要 由 两 个 步骤 来 完成 :(]) 
通过 最 小 化 模型 与 真实 图 像 之 间 的 测度 来 分 解 图 像 , 得 到 
初始 的 因子 与 因子 图 像 -(2) 解 决 模型 中 的 非 唯一 问题 。 


我 们 通过 a 散 度 来 测定 图 像 7 与 模型 CF 之 间 的 偏 
差 ,a 散 度 由 S.Amari 提 出 ,是 诸多 信息 散 度 族 中 的 一 
种 ,用 于 测量 两 个 数据 分 布 p(x) 和 g(x) 之 间 的 偏差 ,其 
定义 为 : 
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D.(p.q)- A *(-2)q,-p/q; "),a eR (4) 


结合 信息 散 度 的 特点 ,本 文 提出 最 小 化 基于 o 散 度 
测度 的 目标 函数 来 估计 因子 与 因子 图 像 C: 
al, +(1 -OCF), 
f.(C,F)= TE 可 y| LCF), IL eR (5) 
基于 a 散 度 的 动态 PET 图 像 的 因子 分 析 模 型 可 表 
示 为 以 下 凸 优化 问题 


(6,5) arg amm fap) (6) 

我 们 采用 一 种 迭代 的 方式 去 解决 上 式 , 具 体 算法 
如 下 : 

i 初始 化 C,F: 设 置 C 中 元 素 为 0-1 之 间 的 随机 数 ， 
再 根据 C 估 计 初 始 ; 

iif: 


(a)C-step: 根 据 当 前 的 动态 图 像 T 和 Fr 估计 CC; 
(b) 将 C 中 的 小 于 零 的 元 素 置 为 0:; 

(c)F-step: 根 据 当前 的 动态 图 像 T 和 CC 估计; 
(d) 将 正中 的 小 于 零 的 元 素 置 为 0:; 

(e) 检 查 是 否 收敛 , 若 不 收敛 , 则 返回 步 又 ai; 

过. 归 一 化 因子 图 像 C, 再 根据 缩放 的 比例 缩放 无; 


Of (C, F P 
其 中 ,在 C-step 中 , 令 HCP -0 得 到 
ik 
T j 
GEN EL o 
, 23 kf | Ta 
Ca = ; 
YF, 
f 
在 F-step 中 , 令 =0, 得 到 
Af 
I, y 
P C, : 
a » fes) 
Fu I 
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1.3 唯一 性 约束 

容易 发 现 , 式 1 描述 的 因子 模型 的 结果 在 数学 上 并 
不 是 唯一 的 呈 ,因子 和 因子 图 像 的 非 唯一 性 可 以 通过 下 
式 说 明 

I-(CR(R'F)-C'F' (7) 

Her, CBLE Xo ECSEBS DSL EUH B [RTL EST b 
[^ at E —^4- np RABIA R, Elie: EEIEIEE 
模型 描述 的 因子 图 像 C 和 已。 因此 ,我 们 必须 找到 合理 
的 约束 条 件 确保 结果 的 唯一 性 。 

A 传统 方法 :最 小 化 因子 图 像 之 间 的 重合 程度 约束 

sitek.A 假 设 因 子 图 像 之 间 的 重 半 程度 最 小 (MSO,， 
minimal structure overlap) 的 作为 唯一 性 约束 条 件 ， 
子 间 重 准 程 度 的 计算 公式 为 

Bn dCR] ICal 


p-lg-ptli-l Y (cR) Y (CR 
il yw rem ic 


(8) 


惩罚 项 f(CR,R 站 保证 因子 和 因子 图 像 的 值 为 非 
负 , 可 用 下 式 表示 : 
f(R)- XX H(CR.) + Y (acp, 


lk-1 


0, x20 o, 
人 [o x<0 

EME H b ER Gur bif) , 求 得 R 的 最 优 解 ,并 获 
得 最 终 的 因子 图 像 (CR) 和 代表 每 个 组 织 活 度 曲 线 (RF) 
的 具有 实际 生理 意义 的 因子 。 

B. 新 方法 :基于 动力 学 聚 类 的 约束 

由 于 受到 PET 图 像 分 辨 率 和 部 分 容积 效应 的 影 
响 , 心 肌 信号 经 常会 摊 杂 着 10%~15% 的 血液 信号 ,最 小 
化 因子 岁 像 之 间 的 重 有 到 程度 约束 未 能 考虑 到 这 一 事 
实 。 因 此 ,本 文 提出 一 种 基于 动力 学 的 时 间 活 度 曲 线 聚 
类 (KB, Kinetic Based) 的 先 验 约束 来 解决 因子 分 析 模 
型 中 的 非 唯一 问题 。 

在 动态 PET 人 研究 中 ,每 一 个 像素 的 时 间 行 为 可 以 
通过 TAC 曲线 来 描述 ,因此 ,我 们 基于 像素 的 动力 学 对 
图 像 进 行 分 类 。 我 们 的 目标 是 通过 动态 序列 图 像 的 
TAC 曲线 的 形状 和 级 数 对 图 像 进 行 分 类 ,使 得 像素 
TAC 在 同一 类 是 同 质 的 ,而 不 同类 之 间 是 异 质 的 。 假 设 
动态 图 像 1 具 有 KK 个 不 同 特征 的 类 ,模糊 C 均 值 聚 类 算 
法 基于 最 小 化 目标 国 数 


J-2Y Fai] -vl (10) 

其 中 | e J| 代表 欧 氏 距离 ,7 表示 图 像 7 中 第 /个 像素 
的 TAC 曲线 ,w 是 第 kt 个 聚 类 的 中 心 的 TAC,mw 表 示 第 / 
个 像素 属于 第 类 的 来 属 度 ,q(g>1) 是 模糊 参数 。 

目标 函数 的 解 可 以 通过 如 下 迭代 式 得 到 | 

i 设 置 初始 的 类 数 尖 和 停止 准则 e; 

计 随 机 初始 化 模糊 隶属 度 矩 阵 D7" 包 含 每 一 个 单独 
uÉ; 

ii D EBENE USE 5-0; 

iv. UAR bv” 

buo P 


o. j 


VREE Un 

O+) c URAN did 

2 =]1 $ 

". SE) 

vi MUR max [U^ - Uc Wr LE; UL UE bb 
LEES iv2E 

将 聚 类 中 心 w 通 过 式 (11) 作 为 像素 动力 学 先 验 信 
息 应 用 到 因子 分 析 唯一 性 约束 当中 ， 


fac) 2] GC Py -vl (11) 


过 最 小 化 目标 函数 (fiict+5j,) 来 解决 因子 分 析 中 
的 非 唯一 问题 , 求 得 R 的 最 优 解 , 最 终 获 得 因子 图 像 
(CR) 和 代表 每 个 组 织 活 度 曲 线 (R"'F) 的 具有 实际 生理 
意义 的 因子 。 


i 


2 实验 

(1) 计 算 机 仿真 数据 

在 计算 机 仿真 实验 中 ,本 文 使 用 基于 如 图 2 所 示 的 
XCAT 数 字体 膜 ,心脏 区 域 由 血 池 ,心肌 组 织 两 个 部 分 
构成 ,在 这 里 将 心脏 区 域 单独 分 离 出 来 。 模 型 的 体 素 点 个 
数 为 64x64x30, 体 素 点 尺寸 大 小 为 3.27 mmx3.27 mmx 
3.27 mm。 我 们 基于 单 房 室 模型 天 模拟 心脏 区 域 生理 代 


AS dH T 


Lic 


谢 的 功能 成 像 过 程 ,使 用 的 放射 性 核 素 为 *Rb, 半 衰 期 
为 1.27 mins 仿真 过 程 中 , 血 输入 函数 以 及 各 组 织 对 应 
的 动力 学 参数 的 测量 值 均 来 自 美国 约翰 霍 普 金 斯 大 学 
PET 医 学 中 心 的 5 位 志愿 者 。 基 于 单 房 室 模型 血 输 
入 函数 以 及 相应 的 动力 学 参数 (设置 动力 学 参数 K= 
1.4822 mL/min/g, k;-0.3159/min 及 血 体 积分 数 为 V; 
0.3829) 生 成 心肌 组 织 的 时 间 活 度 曲 线 , 各 组 织 时 间 活 
度 曲 线 及 空间 分 布 如 图 3 所 示 。 整 个 动态 PET 扫 描 协 
议 设置 为 :12x5 s,6x10 s 共 计 18 个 时 间 帧 。 将 每 一 帧 
血 池 与 心肌 组 织 的 时 间 活 度 赋值 到 其 对 应 组 织 的 空间 
分 布 位 置 上 ,生成 模拟 心脏 区 域 的 动态 PET 图 像 。 

将 动态 图 像 进行 正 向 投影 得 到 随时 间 变 化 的 无 噪 


图 2 (44D XCAT 数 字体 膜 的 正 位 像 ( 左 :男性 , 右 :女性 ) 及 (B) 女 性 心脏 区 域 水 平 切片 
Fig.2 Anterior views of the 4D XCAT digital phantom (A; left: male, right: female) and the female 


horizontal plane of myocardial region (B). 


x10 
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图 3 血 池 与 心肌 组 织 的 时 间 活 度 曲 线 及 空间 分 布 位 置 ( 横 截 面 第 14 层 ) 
Fig.3 Time-activity curve and spatial distribution position of blood pool and 
myocardial tissue (14th slice of the horizontal plane). 


投影 数据 。 投 影 数 据 仿真 过 程 中 产生 的 总 光子 计数 设 
置 为 Sx10'。 在 投影 数据 中 加 入 泊 松 噪声 便 可 得 到 仿真 
的 加 噪 投影 数据 。 本 文采 用 FBP 重 建 算法 对 加 咯 的 投影 


数据 进行 重建 。 得 到 重建 图 像 之 后 ,用 FWHM=4 mm 
的 高 斯 滤波 器 对 重建 出 的 动态 图 像 进 行 去 噪 。 本 研究 
仿真 实验 均 基 于 该 仿真 数据 。 
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(2) 对 照 方法 及 评价 指标 

在 模型 的 求解 中 的 第 一 步 中 ,引入 了 c 散 度 作 为 测 
度量 化 模型 与 真 值 图 像 之 间 的 距离 ,我 们 将 其 与 传统 的 
最 小 二 乘 测 度 进 行 比 较 。 本 文采 用 归 一 化 误差 58 评价 
不 同 测度 在 每 一 帧 的 剩余 信号 s, 用 归 一 化 误差 2 评价 不 
同 测度 评价 所 有 帧 的 剩余 信号 。 


ô= 1 
Ni 
1 T N 
533085 
6,= (12) 
N72 


为 了 验证 本 文中 动力 学 聚 类 在 解决 非 唯一 性 问题 
的 准确 性 ,我 们 将 基于 动力 学 聚 类 约束 的 因子 分 析 模型 
与 传统 的 最 小 化 因子 间 重 蕉 程度 (MSO) 的 模型 进行 比 
较 。 除 了 视觉 评价 外 ,采用 量化 指标 RMSE 对 因子 与 仿 
真 的 组 织 活 度 曲线 分 布 进行 比较 : 


14.0%F 
12.0% 
10.0%| - 

8.0% 


6.0% 


4.0% 


0 2 4 6 8 10 12 14 16 18 
Frame 


图 4 TAWE TREE MARES, 
Fig.4 Normalization error ð of the model in each 
frame under different measures. 


3.2 唯一 性 约束 

得 到 初始 因子 后 ,用 最 小 化 因子 间 重 本 程度 
(MSO) 和 基于 动力 学 聚 类 (KB) 的 两 种 约束 方法 解决 因 
子 分 析 模 型 中 的 非 唯 一 问题 。 图 6、 图 7 分 别 给 出 了 对 
我 也 燥 比 (SNRJ)=10 的 动态 序列 用 两 种 不 同 约束 方法 
(模型 的 测度 为 a 散 度 测度 , 且 a=25) 所 得 到 的 因子 和 
子 图 像 ,由 图 可 知 ,所 得 出 的 因子 1 与 因子 2 分 别 对 应 于 
血 池 TAC 与 心肌 TAC, 所 得 到 的 因子 图 像 1 与 因子 图 像 
2( 图 中 绘制 的 是 第 14 层 横断 面 图 ) 也 对 应 于 心脏 的 血 
和 心肌 组 织 , 对 于 血 池 的 活 度 曲 线 , 两 种 方法 所 得 到 
因子 的 精确 程度 无 显著 差别 ,而 对 于 心肌 活 度 曲 
线 , 用 KB 约束 得 到 的 因子 比 MSO 约束 得 到 的 因子 精 
确 度 高 。 


I LX 
cmo 


Yr, norm(F) - norm(F,,,) | 
RMSE = 一 
YWnorm(F,,) 
其 中 F 代 表 因子 曲线 ,Ps 代表 各 组 织 的 活 度 曲线 。 
Wie URS EST TRITRTISS norm G7 RE F ET TUH— 
化 , 即 norm(F) SF . 


x100% (13) 


J 


3 结果 
3.1 o 散 度 与 最 小 二 乘 测度 建 模 的 精确 性 的 比较 

通过 最 小 化 模型 与 真实 图 像 之 间 的 测度 得 到 初始 
因子 和 因子 图 像 后 ,可 用 剩余 信和 号 来 评价 不 同 测度 的 优 
劣 。 图 4 是 最 小 二 乘 测度 与 x 散 度 测 度 (o=1,10,20,30， 
40,50) 所 得 出 的 模型 归 一 化 误差 6,, 图 5 是 动态 图 像 所 
有 帧 的 归 一 化 误差 5 随 a 的 变化 情况, 由 图 可 知 , 仅 当 o= 
1 时 a 测度 与 最 小 二 乘 测度 的 归 一 化 误差 相等 , 当 a 取 其 
它 值 时 ,a 散 度 测度 的 结果 的 归 一 化 误差 均 优 于 最 小 二 
乘 测度 。 


8.50% f 
— a divergence measure 


8.40% f — Least square measure 
8.30% 上 
8.20% 上 
8.10% 上 
8.00% 上 


7.90% 上 


7.80% 


5 10 15 20 25 30 35 40 45 50 
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图 5 不 同 测度 下 模型 的 归 一 化 误差 8 
Fig. Normalization error à; of the model under 
different measures. 


3.3 o 值 的 选取 

为 了 进一步 验证 a 散 度 在 本 模型 中 的 有 效 性 基 选 取 
最 优 的 a 值 ,本 文选 取 不 同 的 a 值 在 两 种 约束 模型 下 进行 
因子 分 析 , 为 避免 实验 受 噪声 的 随机 性 影响 ,针对 每 种 
方法 运行 30 次 求 平均 值 ,绘制 出 不 同 x 值 在 MSO 和 KB 
两 种 约束 下 各 因子 的 RMSE 图 8, 在 KB 约束 中 ,心肌 因 
子 的 RMSE 基本 保持 稳定 , 当 o=5 时 估计 的 血 池 因子 
RMSE 最 小 ,因此 选取 5 作为 本 文 的 最 优 的 a 值 ;而 在 
MSO 约 束 中 ,两 因子 在 o>12 时 基本 趋 于 稳定 。 选 取 了 
最 优 的 值 后 ,将 a 散 度 测度 与 最 小 二 乘 测度 在 两 种 不 同 
约束 下 作 比 较 见 表 1, 结 果 说 明 a 散 度 的 RMSE 比 最 小 
二 乘 的 低 , 即 该 测度 更 能 提取 出 更 准确 的 结果 。 此 外 ， 
由 图 8 可 知 无 论 c 取 何 值 ,用 KB 约束 模型 得 出 的 两 因子 
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图 6 用 MSO 约束 所 得 因子 及 因子 图 像 
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Fig.6 The factors and factor images obtained by MSO constraints. 
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图 7 用 KB 约束 所 得 因子 及 因子 图 像 


Fig.7 Factors and factor images obtained by KB constraints. 


的 平均 RMSE 曲线 位 于 MSO 的 下 方 , 这 亦 说 明了 动力 
学 聚 类 约束 要 优 于 小 化 因子 间 重 闭 程 度 约束 。 
3.4 模型 对 骂 声 的 敏感 度 

为 了 评价 噪声 对 模型 的 影响 ,我 们 通过 在 图 像 的 投 
影 域 中 加 入 不 同 规模 的 噪声 ,重建 出 不 同 信 噪 比 的 图 
像 。 表 2 表示 对 不 同 信 噪 比 的 动态 序列 ,用 MSO 和 KB 
两 种 不 同 约束 所 得 到 的 各 因子 的 RMSE。 由 表 可 知 , 当 
噪声 逐渐 增 大 时 ,因子 的 RMSE 也 相应 增 大 ,估计 的 血 
池 的 时 间 活 度 曲 线 比 心肌 活 度 曲线 精确 。 


4 讨论 

a 散 度 测度 的 目的 是 将 动态 图 像 做 一 个 初步 的 分 
解 ,由 于 分 解 的 结果 具有 非 唯 一 性 ,所 以 根据 聚 类 中 心 
作为 先 验 信息 将 其 中 一 组 解 线 性 变换 为 具有 实际 生理 


不 能 作为 最 终 的 TAC ,因为 聚 类 与 聚 类 中 心 过 于 粗糙 ， 
不 能 较 精 确 地 反应 真实 的 动态 网 像 ,尤其 是 对 于 噪声 较 
大 的 动态 网 像 。 

在 加 入 唯一 性 约束 时 ,我 们 引入 一 个 惩罚 项 来 保 
证 我 们 得 到 的 结果 是 非 负 的 ,对 于 两 个 正则 化 参数 b, 和 
bs, 实 验 发 现 当 b,,bye[10, 500] 时 ,结果 无 明显 差异 ,鉴于 
此 在 本 实验 中 将 b,b; 均 设置 为 100。 对 于 不 同 的 数据 ， 
我 们 可 以 参考 本 数据 选取 合适 的 bp 值 , 耕 结 果 出 现 负 值 ， 
我 们 只 需要 适当 地 增加 2 的 值 保证 结果 是 非 负 的 即 可 。 

从 该 心脏 体 膜 中 提取 的 结果 我 们 我 发 现 ,在 同样 规 
模 噪 声 的 情况 下 ,估计 出 的 血 池 的 时 间 活 度 曲线 比 心肌 
组 织 的 活 度 曲线 要 精确 ,主要 原因 是 血 池 的 体 素 数 比 心 
肌体 素数 多 , 受 部 分 容积 效应 及 噪声 的 影响 较 小 。 在 心 
肌 壁 的 某 些 部 分 仅 有 2~3 个 体 素 的 厚度 , 极 容易 受到 部 


意义 的 时 间 活 度 曲 线 。 值 得 注意 的 是 ,动力 学 聚 类 结 


分 容积 效应 及 噪声 的 影响 降低 结果 的 精确 度 。 
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—— KB blood factor 
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图 8 不 同 o 值 在 MSO 和 KB 两 种 约束 下 各 因子 的 RMSE 
Fig.8 RMSE of the obtained factors with MSO constraints and KB 
constraints undera divergence measure. 


表 1 两 种 测度 下 应 用 MSO 约束 和 KB 约束 的 所 得 各 因子 的 
RMSE 

Tab.1 RMSE of the obtained factors with MSO constraints and KB 
constraints under two measures 


Least square measure adivergence measure (0=5) 


Item 

MSO KB MSO KB 
Blood 4.0796 1.696 2.4396 0.52% 
Myocardium 4.89% 1.48% 4.43% 1.45% 


表 2 不 同 信 噪 比 下 各 因子 的 RMSE 
Tab.2 RMSE of various factors under different SNRs 


No noise SNR-10 SNR-5 
Item 
MSO KB MSO KB MSO KB 
Blood 0 0 2.90% 0.64% 8.96% 1.3396 
Myocardium 0 0 44499 1.44% 11.32% 0.46795 


Notes: MSO representsMinimal Spatial Overlap constraint method and 
KB represents Kinetic Basedconstraint method. 


本 研究 提出 一 种 基于 动力 学 聚 类 与 o 散 度 测度 的 动 
三 维 心 脏 PET 图 像 因 子 分 析 研 究 模 型 ,实现 了 对 心脏 
PET 图 像 各 组 织 的 时 间 活 度 曲 线 进 行 精确 提取 。 传 统 
的 因子 分 析 模 型 利用 最 小 二 乘 模型 ,用 最 小 化 因子 图 像 
se ed 

过 最 小 化 动态 图 像 与 因子 分 析 模 型 的 wx 散 度 测度 ， 
ae 同 的 c 值 选取 比 最 小 二 乘 适 合 的 测度 得 到 初 
始 因 子 及 因子 图 像 , 并 在 旋转 因子 时 引入 PET 图像 像素 
的 动力 学 聚 类 先 验 信 息 ,我们 将 模型 应 用 于 基于 XCAT 
的 仿真 数据 上 ,通过 实验 分 析 发 现 ,用 a 散 度 测度 以 及 基 
于 动力 学 聚 类 的 先 验 约束 均 对 组 织 活 度 曲 线 的 提取 的 


态 


准确 性 有 积极 的 作用 ,在 高 噪声 下 (SNR=5) 仍 能 准确 的 
提取 出 各 组 织 的 时 间 活 度 曲 线 , 用 视觉 评价 及 量化 评价 
等 方法 优势 明显 。 然 而 ,本 方法 在 旋转 因子 时 引入 正则 
化 参数 无 法 目 适 应 选择 ,在 下 一 步 工 作 中 将 进一步 研究 


正则 化 参数 自 适应 选择 问题 。 
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